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The stacked triangular lattice has the shape of a triangular prism. In spite of being consid- 
ered frequently in solid state physics and materials science, its percolation properties have received 
few attention. We investigate several non-universal percolation properties on this lattice using 
Monte Carlo simulation. We show that the percolation threshold is p^"'''^ = 0.186 02 ± 0.000 02 for 
bonds and pf^^ = 0.262 40 ± 0.000 05 for sites. The number of clusters at the threshold per site 
is n^^^'i = 0.284 58 =b 0.000 05 and nf ^ = 0.039 98 ± 0.000 05. The stacked triangular lattice is a 
convenient choice to study the RGB model [Sci. Rep. 2, 751 (2012)]. We present results on this 
model and its scaling behavior at the percolation threshold. 

PACS numbers: 64.60.ah, 64.60.al, 89.75.Da 



I. INTRODUCTION 

Percolation is a classical model in statistical physics, 
exhibiting a continuous transition between a macroscop- 
ically disconnected and a connected state [1, 2]. The 
elements of a lattice (bonds or sites) are present with 
probability p and absent with probability 1 — p. Sites 
connected by bonds constitute clusters. Above a critical 
value p > Pc^ there is at least one cluster connecting two 
sides of a lattice with free boundary conditions. At the 
critical point p = Pc^ scale- free behavior is observed. 

Many important properties of the percolation transi- 
tion, such as critical exponents, correction-to-scaling ex- 
ponents, and amplitude ratios are universal, and so in- 
dependent on the microscopic lattice details, while oth- 
ers do depend on the considered lattice [1, 3-6]. The 
non-universal properties include in particular the critical 
value Pc of the occupation probability p, i.e., the percola- 
tion threshold of the lattice, and the number of clusters 
per site Uc at the threshold. For applications of percola- 
tion theory, it is important to know these non-universal, 
lattice dependent, properties with precision [2, 7]. 

Here, we investigate the percolation properties of the 
stacked triangular lattice, which has the shape of a trian- 
gular prism (see Fig. 1). Although the stacked triangular 
lattice is used frequently for models in solid state physics, 
for example in studies of antiferromagnets [8-11] or hard- 
core bosons [12, 13], its percolation properties have re- 
ceived far less attention; to our knowledge, only the cal- 
culation of its percolation threshold has been carried out, 
with spanning and wrapping probabilities [14-16], giving 
^bond ^ Q ;Lg5 9 ^ Q QQQ 2 and ^ = 0.262 3 ± 0.000 2. 

The geometry of the stacked triangular lattice is a con- 
venient choice to study the recently introduced RGB 



model [17]. This model is inspired by the challenge 
of sharing reservoirs, for example, when aquifers or 
petroleum reservoirs are concerned [18-21]. The RGB 
model is based on the intuitive fact that partitioning a 
volume into two parts requires one division surface, while 
three parts of a volume are in contact along a simulta- 
neous boundary consisting of lines. It has been shown 
that this boundary is a fractal set consisting of a sin- 
gle thread of fractal dimension 1.55 ± 0.03 spanning the 
medium surrounded by a cloud of loops of fractal dimen- 
sion 1.69 ±0.02 [17]. 

The RGB model generalizes the model of fracturing 
ranked surfaces [22] and the watershed model [23-25] to 
three dimensions and three basins. It is then closely re- 
lated to loopless percolation and the random fuse model 
in strong disorder [22, 26, 27], which exhibit the same 
fractal dimension = 1-2168 ± 0.0005 [28] also ob- 
served for optimum paths, minimum spanning trees [29- 
32], optimum-path cracking [33, 34], and fractal cluster 
boundaries of discontinuous percolation models [35, 36]. 
Conceptually, there is also a relation to the water re- 
tention model, where a fractal depletion zone arises, for 
which, however, a higher fractal dimension of about 5/4 
has been reported [37, 38]. 

The organization of this work is as follows. Section II 
defines the RGB model on the stacked triangular lattice. 
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FIG. 1. (Color online) Sketch of a stacked triangular lattice 
of height L = 6 and lateral size 5 = 3. 
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along with the first results. In Section III, the perco- 
lation thresholds for site and bond percolation on the 
stacked triangular lattice are determined with improved 
precision. These results are used in Section IV to ana- 
lyze the percolation properties and the RGB model at the 
percolation threshold. Conclusions are drawn in Section 
V. 



II. RGB MODEL 

The stacked triangular lattice of height L and lat- 
eral size B consists of L stacked layers of triangular 
lattices of length B (see Fig. 1). Therefore, it has 
Ns,/\ = L{B + B'^)/2 sites and, with free boundary con- 
ditions, Nb,a = B{ABL - B - 2L-l)/2 bonds. 

On this lattice, the RGB model, introduced recently 
in Ref. [17], can be formulated as follows. Initially, all 
Nb^a bonds are unoccupied, such that there are Ns^a 
clusters of unitary size. Some special sites (typically in 
the boundaries) are selected and labeled as red, green, 
or blue, such that every site has at most one color. We 
consider an initial site coloring as shown in Fig. 2(a): 
The three rectangular faces of the prism are divided into 
three parts of the same area as shown in the figure, such 
that sites in the same face have the same color. 

Starting from this setup, there are two equivalent ways 
to simulate the RGB model. Similar to invasion percola- 
tion [39, 40], one can consider that the lattice is invaded 
simultaneously from all initially colored sites, with the 
constraint that clusters of different colors can not merge 
[17]. Here, as an alternative, we consider the following 
procedure, which can be defined in the language of ran- 
dom percolation. Starting from the same initial condi- 
tion, the following steps are iterated until all Nb^/\ have 
been selected: 

(1) Select a bond uniformly at random from the list of 
bonds that have not been considered; 

(2) Consider the colors of the two sites that this bond 
would connect: 

(2.1) If at least one of the two sites has no color or both 
colors are the same, the bond is occupied, merging the 
two sites into the same cluster or closing a loop. In the 
case that one of these clusters has color c while the second 
one is uncolored, all sites in the second cluster become 
colored with c, and the two clusters merge. 

(2.2) Else, i.e., in case occupying the selected bond would 
merge two clusters with different colors, the bond is not 
occupied and it is identified as a bridge bond [22]. 

After all Nb^/\ bonds have been selected, each one is 
either occupied or has been identified as bridge bond, 
and the Ns^a sites in the lattice are split into three clus- 
ters, corresponding to the three colors: red, green, and 
blue [see Fig. 2(b)]. These three clusters are compact, 
but their boundaries are highly corrugated. To quantify 
this observation, we define an RGB site as a site with 
at least three of its neighboring sites carrying distinct 
colors. For a given configuration of colored sites, we can 




FIG. 2. (Color online) (a) Sketch of the initial conditions 
for the RGB model on the stacked triangular lattice. The 
cubes indicate the initial coloring of the sites and the sites not 
located in the rectangular faces of the prism initially have no 
color, (b) Snapshot of a typical final configuration, showing 
a slice parallel to the triangular faces of the lattice (z = 0). 
The RGB sites are marked in white and indicated with the 
white arrow. The lattice size is L = 5 = 128. 



then measure the total number of RGB sites, Mtot, in the 
lattice. In addition, we know the total number of bridge 
bonds, Msur, forming the division surfaces which sepa- 
rate pairs of colored clusters. Finally, as can be observed 
in Fig. 3, due to the special geometry of the stacked tri- 
angular lattice and the chosen initial coloring conditions 
(see Fig. 2), the set of all RGB sites always contains a 
spanning path of RGB sites connecting the two triangu- 
lar faces of the lattice. One can use the burning method 
to determine the length, Mpat, of the shortest path of 
RGB sites connecting the two triangular faces [41]. 

Figure 4 shows the results of measuring Mtot, Mpat, 
and Msur on lattices of different sizes with L = B and 
averaging over many realizations. To keep track of the 
cluster properties, we considered the labeling scheme pro- 
posed by Newman and Ziff [43-45] , related to the Hoshen- 
Kopelman algorithm [46]. The algorithm proposed in 
Ref. [47] was used to generate random numbers. Un- 
less indicated, error bars are smaller than the symbol 
size. For large lattice sizes L the total number of RGB 
sites scales with L as Mtot ^ L^^°^ with fractal dimen- 
sion dtot = 1-69 ± 0.02, as shown in Fig. 4. These results 
are confirmed by the analysis of the local slopes of the 
data and they are compatible with the fractal dimension 
obtained in Ref. [17] using an alternative definition of 
the RGB set. The length of the shortest path of RGB 
sites connecting the two triangular faces of the lattice 
scales as Mpat ~ L"^^""^ with (ipat = 1.54 ±0.03. Using the 
yardstick method [42], we obtain within error bars the 
same result, 1.53 ± 0.05 (see lower inset of Fig. 4). We 
note that dpat is larger than the fractal dimension of the 
minimum path in three-dimensional critical percolation, 
1.3756 ± 0.0006 [48], as well as the one of the optimum 
path in strong disorder 1.44 ± 0.02 [49]. Finally, the size 
of the surfaces separating the three clusters Mgur, mea- 
sured as the total number of bridge bonds, is observed 
for large lattices to behave as M^ur ^ L^^^^. The frac- 
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FIG. 3. (Color online) (a) Visualization of all RGB sites 
of a typical configuration. In blue (dark gray, connecting the 
triangular faces of the lattice) we see the spanning cluster of 
RGB sites, and in green (light gray) the RGB sites not belong- 
ing to the spanning cluster. The lattice size is L = B = 256. 
(b) RGB sites in the shortest path between the two triangular 
faces of the lattice for the same configuration as in (a). 



tal dimension dsur = 2.49 ± 0.03, measured in Fig. 4, 
is compatible with the best known value for the three- 
dimensional watershed fractal dimension, 2.487 ± 0.003 
[28, 50]. 



III. DETERMINING THE PERCOLATION 
THRESHOLD 

The RGB model is known to exhibit a crossover point 
described by a negative exponent at the percolation 
threshold Pc for the lattice [17]. To investigate this phe- 
nomenon, it is necessary to know Pc with good accuracy. 
In the following, we improve the estimates of the thresh- 
olds, obtaining = 0.186 02 ± 0.000 02 for bond and 
pf"" = 0.262 40 ± 0.000 05 for site percolation, compat- 
ible, although with more precision, with previous esti- 
mates [14-16]. 



1.12 




Lattice size L 

FIG. 4. (Color online) RGB sites in the stacked triangular 
lattice. Size of division surfaces Mgur (D), total number of 
RGB sites Mtot (O), and length Mpat (A) of the shortest 
spanning path as a function of the lattice length L — B. 
The solid lines are guides to the eye with slopes 2.49, 1.69, 
and 1.54, respectively. The upper inset shows the rescaled 
number of RGB sites Mtot / L^^^ as a function of the lattice size 
L. Extrapolating the local slopes dL = log2[M(L)/M(L/2)] 
to 1/L ^ 0, we obtain the following estimates for the fractal 
dimensions: dsur = 2.49 ±0.03, dtot = 1.69 ±0.02, and dpat = 
1.54 ± 0.03. Results have been averaged over 10^ samples for 
the smallest system size (L = 5 = 8) and 4.9 x 10^ samples 
for the largest one [L — B — 724). The lower inset shows the 
results of the yardstick method [42] to measure the fractal 
dimension dpat of the shortest path in the spanning RGB site 
cluster. The number n of sticks of length / needed to follow 
the path is plotted as function of / for different lattice lengths 
L [362 (V), 512 (o), and 724 (*)]. The solid line is a guide 
to the eye with slope —1.53. From the yardstick method we 
estimate dpat = 1.53 ± 0.05, which is in agreement with the 
result dpat = 1.54 ± 0.03 obtained in the main plot. Results 
have been averaged over at least 10"^ samples. 



As discussed in Section II, the stacked triangular lat- 
tice is convenient for simulating the RGB model. The 
percolation threshold only depends on the dimension and 
topology of the lattice, but not on the boundary con- 
ditions. Thus, for numerical convenience, to determine 
the percolation threshold and other non-universal prop- 
erties, we consider two joined stacked triangular lattices 
(see Fig. 5), resulting in one simple-cubic lattice with 
one additional bond per square in the x?/-plane. Conse- 
quently, a site with coordinates (xo,^o,^o) is connected 
not only to its six nearest neighbors, as in the cubic 
lattice, but also to the sites at (xq + 1, + 1, ^o) and 
(xo — 1, ^0 — 1, ^o)- This lattice, with periodic boundary 
conditions, has Nb = 4L^ bonds and Ns = sites. 

p denotes the fraction of occupied bonds or sites. We 
analyze the lattice size dependence of several threshold 
estimators (see Fig. 6 and Table I and II). In a lattice 
of Ns sites, the order parameter Pqo of the percolation 
transition is defined as the fraction of sites in the largest 
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FIG. 5. (Color online) Sketch of a lattice obtained by 
merging two stacked triangular lattices with B — L — 3. The 
bonds of one of the lattices are shown as green (bold) lines. 
The resulting topology is a simple-cubic lattice with one addi- 
tional diagonal bond per square in the x^-plane, with L = 3. 
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FIG. 6. (Color online) Size dependence of the threshold 
estimators pc,j (□), Pc,m (O)^ Pc,x (A), and pc,a (v), for 
bond percolation (main plot) and site percolation (inset). We 
use 1/z/ = 1.145, following Ref. [62]. Results are averages over 
at least 1.5 x 10^ samples. 



TABLE I. Bond percolation threshold estimates for the 
stacked triangular lattice. 



Estimator 


Estimate 


Error 


Pc,J 


0.186 015 


0.000 010 


Pc,M 


0.186 016 


0.000 010 


Pc,x 


0.186 015 


0.000 015 


Pc,z 


0.186 017 


0.000 020 


Pc,a 


0.186 015 


0.000 015 


Pc,i 


0.186 


0.000 1 


Combined 


0.186 02 


0.000 02 



and we use the average position Pc,M of the peak in M2 
for finite lattices as the second estimator. Finally, we also 
consider the probability that at least one cluster wraps 
around the lattice (with periodic boundary conditions) 
[43-45, 63, 64]. We define W^c as the probability that 
at least one cluster wraps the lattice in x-direction and 
Wz as the probability that at least one cluster wraps the 
lattice in z-direction. Since for the stacked triangular 
lattice, the z-direction is not equivalent to the x- and the 
^/-directions (see Fig. 1 and 5), one in general can not 
expect Wx and Wz to have the same quantitative behav- 
ior. To estimate the percolation threshold, we consider 
the fractions of occupied bonds Pc^x^ Pc,z^ and Pca, where 
the given sample first wraps in x-, z-, or any direction, 
respectively. For large lattices, the wrapping probability 
curves are expected to intersect at values of p close to 
the percolation threshold Pc (see Fig. 7). Therefore, we 
also consider the fraction of occupied bonds Pc,i, where 
two wrapping probability curves for subsequent lattice 
sizes, Wx{p^L) and Wxip^'^L) intersect (see Table III). 
To keep track of wrapping clusters, we followed the pro- 
cedure outlined, e.g., in Refs. [44, 45, 63, 64]. Our final 
estimates of the percolation thresholds of the stacked tri- 
angular lattice are 

^bond ^ Q ;^gg Q2 ^ Q 000 02 (4) 

for bond and 



cluster (of size Smax), 

Poo=Sme.jNs . (1) 

In the thermodynamic limit, at the percolation thresh- 
old, the rate of change of the order parameter diverges 
[51, 52]. Therefore, we consider as one estimator the 
average fraction of occupied bonds or sites Pc, J at which 
the largest change in P^o occurs [36, 53-61]. The second 
moment of the cluster size distribution is defined as, 

M. = ^E4, (2) 

k 

where the sum runs over all clusters and is the number 
of sites in cluster k. The second moment excluding the 
contribution of the largest cluster is, 

M^=M2-sl,jNs , (3) 



pf^ = 0.262 40 ± 0.000 05 (5) 

for site percolation (see Table I and II). 

While recording the threshold estimators Pc^j and 
Pc,M, i-e., the average value of p where the largest change 



TABLE II. Site percolation threshold estimates for the 
stacked triangular lattice. 



Estimator 


Estimate 


Error 


Pc,J 


0.262 43 


0.000 08 


Pc,M 


0.262 42 


0.000 08 


Pc,x 


0.262 42 


0.000 08 


Pc,z 


0.262 4 


0.000 1 


Pc,a 


0.262 40 


0.000 05 


Pc,i 


0.262 4 


0.000 1 


Combined 


0.262 40 


0.000 05 
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FIG. 7. (Color online) Wrapping probabilities for site per- 
colation as functions of the fraction of occupied sites p for 
different lattice sizes L [from right to left: L = 32 (solid red 
line), 64 (long-dashed green line), 128 (dashed blue line), and 
256 (dotted yellow line)]. Main plot: probability Wz that 
at least one cluster wraps the lattice in z-direction. Inset: 
probability Wx that at least one cluster wraps the lattice in 
x-direction. Results are averages over at least 10^ samples. 



TABLE III. Intersection points pc,i of the wrapping proba- 
bilities Wx for bond percolation (a) and site percolation (b). 

(a) 



L 


2L 




Error 


8 


16 


0.185 8 


0.000 5 


16 


32 


0.185 92 


0.000 10 


32 


64 


0.185 99 


0.000 10 


64 


128 


0.186 01 


0.000 10 


(b) 


L 


2L 




Error 


8 


16 


0.261 4 


0.000 5 


16 


32 


0.262 10 


0.000 10 


32 


64 


0.262 33 


0.000 10 


64 


128 


0.262 40 


0.000 10 



in Poo and M2 occurs, we also measured the average val- 
ues of the largest change in Pooip)^ denoted by J, and 
the maximum of M2(p), denoted by M2 (see Fig. 8). 
One observes that for large lattices, the measured quan- 
tities scale as J - L-^-^siO-Oi ^/ _ ^2.04±o.03^ As- 
suming that J - L-^/^ [54] and - L^/^, where 
/3, 7, and V are the critical exponents related with the 
order parameter, the susceptibility, and the correlation 
length, respectively, this is consistent with the known 
values of these ratios for percolation in three dimen- 
sions, p/y = 0.4774 ± 0.0001 and j/u = 2.0452 ± 0.0001 
[7, 62, 65-67]. 
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FIG. 8. (Color online) Maximum change J (□) of Poo and 
maximum of the second moment M2 (O) functions of the 
lattice size L. Main plot: J/L'^-^^^ and M^/L^ ^^^ as func- 
tions of L. Inset: J and M2 as functions of L. Assuming 
J ~ L~^^^ and M2 L^^^, one sees that the exponents for J 
and the maximum of the second moment, ^ jv — 0.48 =b 0.01 
and jy — 2.04 ± 0.03, are consistent with the reported val- 
ues Pjv = 0.4774 ± 0.0001 and 7/zy = 2.0452 ± 0.0001 
[7, 62, 66, 67]. The data is shown for bond percolation. The 
solid lines are guides to the eye. Results are averages over at 
least 1.5 X 10^ samples. 



IV. PERCOLATION AND RGB SITES 

Having determined the percolation threshold we 
use it to measure the number of clusters per site at 
^bond rf^^^^ which for critical percolation is known 
to be a lattice-dependent constant [68-70] and which, in 
addition, can be considered to determine the percolation 
threshold for loopless bond percolation [26, 56, 71]. 

Figure 9 shows the data for nj?°"^ and n^^*®. Extrap- 
olating to the thermodynamic limit, we obtain for bond 
percolation 



^bond 



and 



0.284 58 ± 0.000 05 



0.039 98 ± 0.000 05 



(6) 



(7) 



Tbond 



for site percolation. We remark that n^"''" is 
slightly higher than the simple-cubic lattice value of 
0.272 931 ± 0.000 000 5 [72]. In loopless bond perco- 
lation, bonds that would close a loop, i. e. that would 
connect two sites which are already part of the same 
cluster, remain unoccupied. Therefore, in comparison 
with classical bond percolation, the value of the perco- 
lation threshold differs in general, p|popiess, bond < ^bond^ 
but for every given configuration the obtained clusters are 
identical to the ones in classical percolation. It follows 
that the number of clusters at the threshold is identical, 
^loopless, bond _ ^bond^ Consider a loopless bond perco- 
lation process, starting from an initial condition where 
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FIG. 9. (Color online) Number of clusters at pc per 
site ric as function of the inverse lattice size for bond 
(main plot, □) and site percolation (inset, Q)- Extrapolat- 
ing to the thermodynamic limit, we obtain for bond per- 
colation = 0.284 58 =b 0.000 05 and for the site case 
n^^*^ = 0.039 98 =b 0.000 05. Results are averages over at least 
10^ samples. 



all Nb bonds are unoccupied, p = 0, such that there are 
Ns clusters of unit size. With increasing every bond 
that becomes occupied merges two distinct clusters into 
one tree, since it can not close a loop. If there are pNs 
occupied bonds, the number of clusters is Ns — pNs- At 
the threshold, p = p|popiess, bond^ ^^^^^ 

_ ploopless, bond^^ ^ ^loopless, bond^^ j^gj 



clusters. Therefore, 

loopless, bond 
Pc 



(1 - n^^°''^)Ns/NE 



(9) 



For example, on the two-dimensional square lattice, this 
yields pioopiess, bond _ (7_3y^)/4 0.450962 [71]. 
While on the square lattice Ns/Nb = 1/2, on the stacked 
triangular lattice Ns/Nb = 1/4, as can be seen either by 
considering the modified simple-cubic lattice with Ns ad- 
ditional bonds (see Fig. 5), or by taking the large B = L 
limit of Ns^a/J^b,a- Therefore, we can derive from our 
result in Eq. (6), for the loopless bond percolation thresh- 
old: 



loopless, bond ^ Q^^yg gg ^ q^qqq q2 



(10) 



In a similar language as the one used for the percola- 
tion model, we can understand the evolution of the set 
of RGB sites in the RGB model, by analyzing its be- 
havior with an increasing fraction of selected bonds [17], 
starting from the initial condition shown in Fig. 2(a). To 
analyze the evolution of the RGB model similarly to a 
bond percolation process, we consider the following: For 
bond percolation, p is the fraction of occupied bonds. 
Let us define for the RGB model that p is the fraction 
of bonds that are either occupied or have already been 
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FIG. 10. (Color online) Rescaled total number of RGB sites 
MtotL^ plotted as function of the scaling variable (p — pc)L^ 
with 6* = 1.33 and pc = Pc°'''^ = 0.18602 for different lattice 
sizes L = B [L = 128 (□), 256 (Q), 512 (A), and 650 
(v)]. The RGB site growth exponent is determined to be 
Ct = 1.8 ± 0.3. The solid line is a guide to the eye with slope 
1.8. Results have been averaged over at least 30 samples. The 
inset shows the rescaled number of RGB sites Mtot / L^*°* as a 
function of the fraction of selected bonds p, with dtot = 1.69 
for different lattice sizes L = B. Results have been averaged 
over at least 300 samples. 



labeled as bridge bonds. Then the description of the 
RGB model in Section II starts with p = and the re- 
sults presented there are for the case p = 1. The inset 
of Fig. 10 shows the rescaled total number of RGB sites 
Mtot/^^*°* as function of the fraction of selected bonds p. 
One observes that below the bond percolation threshold, 
Mtot/^^*°* vanishes for large lattice sizes L. Above 
Mtot/^^*°* is finite and the data points for different lat- 
tice sizes overlap, showing that the fractal dimension of 
the set of RGB sites is c^tot = 1-69 ± 0.02. With the dis- 
tance to Pc^ Mtot is found to scale as Mtot ^ {p — Pc)^^ i 
where we obtain C,t = 1.8 ± 0.3 for the RGB site growth 
exponent. 

This peculiar behavior raises the question of how the 
set of RGB sites behaves at the percolation threshold. In 
Ref. [17], it was found that at p = Mtot ^ ^ with 



t = 3 - 2dsur{Pc) 



(11) 



where (isur(Pc) is the fractal dimension of the surfaces 
dividing pairs of colors ai p = Pc. Since the division 
surfaces at the percolation threshold are formed out of 
bridge bonds that would, once occupied, form a span- 
ning cluster, dsur{Pc) is the critical bridge bond fractal 
dimension. Coniglio has shown that this fractal dimen- 
sion equals where u is the correlation length critical 
exponent of percolation [73, 74]. For three-dimensional 
percolation, 1/u = 1.1450 ± 0.0007 [62] and one obtains 
t ~ 0.71. These results imply that while the size of the 
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division surfaces M^ur diverges at the threshold, the num- 
ber of RGB sites Mtot decreases as a power law in the 
lattice size. Thus, in the thermodynamic limit the set of 
RGB sites is empty. These findings can be summarized 
in the following crossover scaling for the number of RGB 
sites: 



Mtot{p,L) = L-'Ftot[{p-Pc)L'] , (12) 



where the scaling behavior of Mtot {p^ L) in p and L im- 
plies = (dtot +^)/Ct [17]. The scaling function Ftot[^] 
fulfills Ftot[^] ^ x^^ for X > 0. In the main plot of 
Fig. 10, the Ansatz in Eq. (12) is confirmed. The quality 
of the scaling is seen from the increase of the overlap re- 
gion of curves for different lattice sizes, with the lattice 
size. 



V. FINAL REMARKS 

Concluding, we have used the special geometry of the 
stacked triangular lattice to study the recently intro- 
duced RGB model [17]. We also calculated more pre- 
cise values for the bond and site percolation thresholds 
of this lattice, which could be useful for models in solid 
state physics using the same lattice. In addition, we de- 
termined the number of clusters per site for critical perco- 
lation on the stacked triangular lattice and investigated 
the scale-free behavior of the RGB model at the criti- 
cal point. In the future, it would also be interesting to 
extend the RGB model to higher numbers of colors and 
dimensions. 
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